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We analyze which momentum component of the gluon field induces spontaneous chiral symmetry 
breaking in lattice QCD. After removing the high-momentum or low- momentum component of the 
gluon field, we calculate the chiral condensate and observe the roles of these momentum components. 
The chiral condensate is found to be drastically reduced by removing the zero-momentum gluon. 
The reduction is about 40% of the total in our calculation condition. The nonzero-momentum 
infrared gluon also has a sizable contribution to the chiral condensate. From the Banks-Casher 
relation, this result reflects the nontrivial relation between the infrared gluon and the zero-mode 
quark. 
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I. INTRODUCTION 

Spontaneous symmetry breaking is one of the most 
significant and universal mechanisms in physics [l]-BJ. 
Although the QCD Lagrangian possesses chiral symme- 
try in the chiral limit, SU(A r /)^x SU (Nf)ft symmetry is 
spontaneously broken into its subgroup SU(iV/)y. Spon- 
taneous chiral symmetry breaking is remarkably impor- 
tant in hadron physics [5- 7]. Also, it is one of the domi- 
nant origins of mass in our world. 

Chiral symmetry itself is symmetry of quarks, not glu- 
ons. However, spontaneous breaking is dynamically in- 
duced by the nonperturbative interaction of gluons. The 
gluon dynamics is inseparably linked with chiral symme- 
try breaking. Our goal is to determine what momentum 
component of the gluon field induces spontaneous chiral 
symmetry breaking. The relation between the eigenmode 
of quarks and chiral symmetry breaking is known as the 
Banks-Casher relation Q. On the other hand, the rela- 
tion between the momentum component of gluons and 
chiral symmetry breaking is nontrivial. It is easy to ex- 
pect the importance of the low-momentum gluon, but 
difficult to predict the detailed relation due to the non- 
perturbative dynamics of the low-momentum gluon. We 
would like to clarify such a relation nonperturbatively 
by lattice QCD. In other words, we quantitatively inves- 
tigate the relevant gluonic energy scale of spontaneous 
chiral symmetry breaking from lattice QCD. 

To analyze the relevant gluonic energy scale, we con- 
sider momentum space of the gluon field. The gluon 
field is described by the link variable in lattice QCD. 
By manipulating the link variable in momentum space, 
we directly analyze how the momentum component of 
the gluon field affects chiral symmetry breaking. The 
obtained energy scale would also be interesting from the 
viewpoint of a connection to other QCD phenomena, such 
as confinement (9j-[T^|. In lattice QCD, the connection 
between confinement and chiral symmetry breaking is 
investigated in the context of phase transition at finite 
temperature [l3l - tl6| . Our analysis is a different approach 
to reveal this connection. 



We calculate the chiral condensate (qq) in lattice QCD. 
The chiral condensate is an order parameter of chiral 
symmetry breaking in the chiral limit. It is nonzero in 
the symmetry-broken phase and zero in the symmetry- 
restored phase. We denote the flavor-averaged chiral con- 
densate in the lattice unit as 



(1) 



where a is the lattice spacing and S q is the quark propa- 
gator. When the quark mass is finite, the chiral conden- 
sate includes the effect of explicit breaking by the quark 
mass as well as spontaneous breaking. To extract the 
chiral limit in lattice QCD, one calculates with several 
quark masses and extrapolates to the chiral limit. 

In this paper, we calculate the chiral condensate in 
SU(3) C quenched and full lattice QCD, and analyze the 
relevant gluonic energy scale of spontaneous symmetry 
breaking. This paper is organized as follows. In Sec. II, 
we explain how to analyze the relevant gluonic energy 
scale in lattice QCD. In Sec. Ill, we show the simulation 
setup of the lattice QCD calculation. In Sec. IV, we 
present the numerical result of the chiral condensate and 
analyze how the chiral condensate is affected by removing 
the high-momentum or low-momentum gluon. Finally, 
Sec. V is devoted to a conclusion. 



II. FORMALISM 

The lattice framework to determine the relevant glu- 
onic energy scale was proposed in Ref. [Trj ] . In this frame- 
work, after artificially removing some momentum compo- 
nent of link variables, one calculates a physical quantity 
and observes the role of the removed momentum compo- 
nent. In doing so, one can determine whether the mo- 
mentum component is relevant or not for the quantity. 
To be self-contained, we briefly introduce the procedure 
in the following. 

Step 1. The SU(3) C link variable U^(x) is generated 
by Monte Carlo simulation. As explained below, the link 
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FIG. 1: The schematic figure of momentum space. The 
shaded regions are the cut regions by the ultraviolet cutoff 
Auv and the infrared cutoff Air. The momentum-space lat- 
tice spacing is a p — 2n/ La. 



variable must be fixed with a certain gauge. In this paper, 
we use the Landau gauge for the numerical calculation. 
In the Landau gauge, the gauge fluctuation is minimized 
and the connection between the link variable and the 
gauge field is straightforward. 

Step 2. The momentum-space link variable E/^(p) is 
obtained by the Fourier transformation, as 



J V si to 



(2) 



where iVsite is the total number of lattice sites. 

Step 3. Some component of U^{p) is removed by in- 
troducing a momentum cutoff. In the cut region, the 
momentum-space link variable is replaced by the free- 
field link variable 
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For example, in the case of the ultraviolet cutoff Atjv, 
the momentum-space link variable is replaced as 



(Vp 2 > A uv ). 
In the case of the infrared cutoff Air , it is replaced as 



(4) 



Step 4. The coordinate-space link variable with the 
momentum cutoff is obtained by the inverse Fourier 
transformation as 



[/»=5»)exp(- 



■jy^yXy). 



(6) 



Since U'{x) is not an SU(3) C matrix in general, U'(x) 
must be projected onto an SU(3) C element U A {x). The 
projection is realized by maximizing the quantity 



ReTr[{C/>)}tE/'(z)]. 



(7) 



Step 5. The expectation value of an operator O is 
computed by using this link variable U A (x) instead of 
Un(x), i.e., {0[U A ]) instead of (O [U]). 

Repeating these five steps with various values of the 
momentum cutoff, we observe the dependence on the 
momentum cutoff. Then, we can determine what mo- 
mentum component of the gluon field is relevant for the 
physical quantity. The framework is applicable to both 
quenched and full QCD in the same way. 

Indeed, this framework is powerful in determining the 
releva nt g luonic energy scale of confinement in quenched 
QCD [13, HI]- By applying this framework to the cal- 
culation of the Wilson loop, it was found that the string 
tension is generated by the infrared gluon below about 1.5 
GeV. By picking up this relevant momentum component, 
the quark-antiquark potential is clearly decomposed into 
the confinement potential and the perturbative potential. 
Hence, the relevant gluonic energy scale of confinement 



was determined to be 



< 1.5 GeV. 



We comment on two points of the framework. The first 
is the gauge fixing in Step 1. In general, since the gauge 
transformation is nonlocal in momentum space, the mo- 
mentum region of the gauge field is a gauge-dependent 
concept. Then, our result would depend on the gauge 
choice. We show the Landau-gauge results in this pa- 
per. Note, however, that one can analyze other gauges 
and the gauge dependence since the framework itself does 
not depend on the gauge choice [13, [H| . 

The second is the projection in Step 4. Although such 
a projection is often used in SU(3) C lattice QCD as a 
workable method, the projection could in principle con- 
taminate the original condition on the momentum cut- 
off. To evaluate how the projection changes link vari- 
ables, we calculate U AA (x) by adopting Steps 2-4 once 
again to U A (x), and check the overlap between them, 
lReTr[{U A {x)}^U AA (x)]. The overlap is found to be al- 
most unity. For example, the deviation from unity is 
about 0.1% at Air = 1.5 GeV. Then, we can expect that 
the projection does not significantly change link vari- 
ables. In fact, we have already reached a steady state 
configuration with the single procedure. 



Ufa) = 



$po (vV < Am) 
%{p) Am). 



(5) 



The schematic figure is shown in Fig. [TJ 



III. SIMULATION SETUP 

The lattice QCD simulations are performed in SU(3) C 
quenched and full QCD. The parameters of gauge con- 
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figurations are summarized in Table Q] For the full QCD 
calculation, we use the dynamical configuration which in- 
cludes the two-flavor staggered quark in NERSC archive 
[IeH . The momentum-space lattice spacing a p is given by 
a p = 2tt/ La, where L is the number of lattice sites in the 
spatial direction. 

To compute the chiral condensate, we adopt the stag- 
gered fermion action, which preserves the U(l) subgroup 
of the full chiral symmetry in the chiral limit. In full 
QCD, we use a single mass for the valence and sea quarks, 
ma = m sca a = 0.010. The corresponding pion mass is 
about 500 MeV and the flavor-averaged chiral conden- 
sate is about (540 MeV) 3 . In quenched QCD, we use the 
quark masses ma = 0.010, 0.015, and 0.025 to extrapo- 
late the chiral limit. 
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IV. LATTICE QCD RESULT 



A. Chiral condensate with the UV cutoff 



FIG. 2: The chiral condensate E = ~a 3 (qq)/Nj with the UV 
cutoff Auv The quark mass is ma — 0.01. The "renormalized" 
chiral condensate is multiplied by the renormalization factor 
Zs. 



First, we show the chiral condensate with the ultra- 
violet (UV) cutoff Auv in Fig. [5J Since there is no sig- 
nificant difference between the quenched and full QCD 
results, we plot only the full QCD result. The right-side 
point at Auv — 12.5 GeV is the result of original lattice 
QCD without the momentum cutoff. 

Although spontaneous chiral symmetry breaking is ex- 
pected to be caused by nonperturbative gluons, the chiral 
condensate is drastically changed by the UV cutoff. How- 
ever, as shown below, this is mainly because the chiral 
condensate is a renormalization-group variant and UV- 
diverging quantity. It is dressed by perturbative gluons 
and its value strongly depends on the UV regularization. 
In standard lattice QCD, the perturbative contribution 
is several orders of magnitude larg er than the nonpertur- 
bative core of the condensate 



where 



20] 



To estimate the effect of renormalization, we calculate 
a renormalization factor, so-called a Z-factor, nonpertur- 
batively (2j [22j. The renormalization factor Zo(k) is 
determined from the amputated Green function of the 
quark bilinear operator O. The renormalization condi- 
tion is imposed as 



z {k)z- l {k)T {k) = i, 



(8) 



TABLE I: The parameters of full and quenched lattice QCD 
configurations. The dynamical quark mass m sca , the con- 
figuration number N con {, the lattice spacing a, and the 
momentum-space lattice spacing a p are listed. 







Volume 


m sca a V C onf 


a [fm] 


a p [GeV] 


Full 


5.7 


16 3 x 32 


0.01 24 - 49 


0.098 


0.79 


Quenched 


6.0 


32 4 


10 


0.100 


0.39 



T (k) = 

G (k) = 
S q (k) = 



1 



tr[^ 1 (fe)Go(fc)5- 1 (fc)P ( 



16N L 
(q(k)Oq(k)) 
(q(k)q(k)). 



(9) 
(10) 
(11) 



Po is the appropriate projection operator. The wave- 
function renormalization factor z\ (k) of the quark 
field is obtained from the conserved vector current, i.e., 
Zv(k) — 1. Note that k is the momentum of the quark 
field, not the momentum of the gluon field. 

We calculate the renormalization factor Zs(k) of the 
scalar operator, and plot the renormalized chiral conden- 
sate Zg(5 GeV) x E in Fig. [2] The renormalized chiral 
condensate is almost independent of the UV cutoff. As 
the UV gluon is removed by the UV cutoff, the bare chi- 
ral condensate approaches the renormalized one. This 
means that the drastic change by the UV cutoff is well 
explained in terms of renormalization. 



B. Chiral condensate with the IR cutoff 

Second, we analyze the chiral condensate with the in- 
frared (IR) cutoff Ajr. We show the full QCD result in 
Fig.[3]and the quenched QCD result in Fig. 2) The quark 
mass is ma = 0.01 in both calculations. In the case of 
the IR cutoff, the chiral condensate does not show the 
drastic change corresponding to renormalization. Then, 
we expect the physical contribution to spontaneous chiral 
symmetry breaking instead of an artifact of renormaliza- 
tion. 

When the IR gluon is removed, the effective quark 
mass would be reduced especially at a large distance. 
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FIG. 3: The full QCD result of the chiral condensate E = 
—a 3 (qq)/Nf with the IR cutoff Air,. The lattice volume is 
f6 3 x 32, and the quark mass is ma = 0.01. PBC and APBC 
mean periodic and antiperiodic boundary conditions, respec- 
tively. 
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FIG. 4: The quenched QCD result of the chiral condensate 
with the IR cutoff. The lattice volume is 32 4 , and the quark 
mass is ma — 0.01. The notation is the same as Fig. [3] 



Thus, we must pay attention to the finite- volume ef- 
fect in Air > 0, even though our lattice volume is large 
enough at Air = 0. We estimate the finite- volume effect 
by changing boundary conditions of the quark propaga- 
tor [H]. In Fig. [3] and Fig. H PBC and APBC mean 
periodic and antiperiodic boundary conditions, respec- 
tively. Since the result is independent of the boundary 
conditions if the lattice volume is large enough, the dif- 
ference between these data should be understood as the 
finite-volume effect. As seen from Fig. [3l the 16 3 x 32 
lattice of full QCD suffers from the finite-volume effect 
in Air > 1.0 GeV. From Fig. 21 the finite- volume effect 
is fairly small for the 32 4 lattice of quenched QCD, al- 



though it gradually grows in Air > 1.5 GeV. 

Both in Fig. and Fig. |H the chiral condensate sud- 
denly gets small around Air = 0. This jump around 
Air = is caused by cutting only the zero-momentum 
link variable £^(0). Despite the change at a single point 
p 2 = 0, the chiral condensate is about 40% reduced. Such 
a large change is not observed in removing other low- 
momentum components. Therefore, the zero-momentum 
gluon is special and it possesses a major contribution 
to the chiral condensate. Note that "zero momentum" 
on momentum-space lattice corresponds to the deep- 
in the con- 



< a, 



infrared region which is roughly i 
tinuum. 

In large Air, since the lattice volume of full QCD 
is not large enough, we analyze the quenched QCD re- 
sult in Fig. [4] When the "nonzero-momentum" gluon of 
\pp^ > a p is removed by the IR cutoff, the chiral con- 
densate gradually decreases. Thus, not only the "zero- 
momentum" gluon but also the "nonzero-momentum" 
gluon contributes to the chiral condensate. The chiral 
condensate continues to decrease even in Air > 1.5 GeV. 
Although it is difficult to perform an accurate analysis in 
large Air due to the finite-volume effect, we can see that 
the chiral condensate is also affected by the gluon in the 
intermediate-momentum region. 



C. Chiral extrapolation 

Next, we consider the chiral extrapolation of the chiral 
condensate. When the bare quark mass m is small, the 
chiral condensate is expanded as a function of m, as 



S(m) = S(0) + maE'(0) 



(12) 



where S'(m) = dY,(m)/dma. E(0) represents sponta- 
neous chiral symmetry breaking in the chiral limit. We 
fit the quenched QCD result by the linear extrapolation 
function E(0) + maS'(0). The fitting result is shown in 
Fig. [5]and Table HI Note that the data of "Air - 0.1 
GeV" corresponds to the smallest IR cutoff, which cuts 
only the zero-momentum link variable, and so the value 
"0.1 GeV" itself is not so meaningful. 

As stated above, when the zero-momentum gluon field 
is removed, the chiral condensate is largely changed. S(0) 
is about 40% reduced and S'(0) is about 30% reduced. 
As for the nonzero-momentum gluon, the extrapolating 
line moves down parallel by the infrared cutoff. E(0) is 
gradually reduced and E'(0) is almost unchanged. This 
indicates that the nonzero-momentum gluon has small 
but finite contribution to spontaneous chiral symmetry 
breaking. 

In Fig. [5l our result suggests another interesting pos- 
sibility. At least within the present numerical accuracy, 
the chiral condensate in the chiral limit remains finite at 
Air = 1.5 GeV, which is the relevant gluonic energy scale 
of confinement. If this is true, this means that the gluonic 
energy scale of spontaneous chiral symmetry breaking is 
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FIG. 5: The chiral extrapolation of the chiral condensate £ = 
—a 3 (qq)/Nf. Air ~ 0.1 GeV corresponds to the cutoff for the 
zero-momentum link variable. 



larger than that of color confinement at zero tempera- 
ture. Unfortunately, however, we cannot make a decisive 
statement due to systematic error of the chiral extrap- 
olation. For more conclusive answer, we need the full 
QCD calculation very close to the chiral limit, while the 
finite-volume effect is severely crucial in large Air and 
small to. 



V. CONCLUSION 

We have analyzed which momentum component of the 
gluon field induces spontaneous chiral symmetry break- 
ing. The most dominant contribution is given by the 
"zero-momentum" gluon, which roughly corresponds to 
the deep-infrared region of \pp* < a p in the contin- 
uum. Not only zero-momentum but also the "nonzero- 
momentum" gluon of \fp^ > a p possesses a sizable con- 
tribution. While we cannot precisely determine the upper 
limit of the relevant momentum component of the gluon 
field due to the finite- volume effect, its relevant momen- 
tum component seems to be broadly distributed to the 



TABLE II: The fitting result of the chiral extrapolation in 
Fig. [5] The extrapolation function is E(0) + ma£'(0). 

Am S(0) X'W 

0.00639(81) 1.269(53) 

- 0.1 GeV 0.00380(27) 0.933(16) 

1.5 GeV 0.00200(7) 0.948(4) 

1.8 GeV 0.00155(2) 0.929(1) 



intermediate-momentum region. 

The zero-momentum gauge field corresponds to a 
spatially-uniform gauge background. In general, the non- 
Abelian gauge field could have a nontrivial effect even 
in spatially-uniform case, unlike the Abelian gauge field. 
Our result actually suggests that the zero-momentum 
gauge field contributes to the chiral condensate. Note, 
however, that it is nontrivial whether spontaneous chiral 
symmetry breaking occurs only by the spatially-uniform 
gauge background. 

The Banks-Casher relation states that the chiral con- 
densate is related to the spectral density p(X) of the Dirac 
operator as 

(qq) = -7T P (0) (13) 

in the chiral limit Q . The spectral density of the Dirac 
operator is given in infinite volume as 



p(X)= lim -Y,S(\-\ k ), (14) 

k 

and the eigenvalue of the Dirac operator is i\k- The zero 
mode of quarks is directly related to spontaneous chiral 
symmetry breaking from this relation. In contrast, the 
gluon field is nontrivially related to spontaneous chiral 
symmetry breaking. Our result presents the connection 
between the momentum component of gluons and the 
zero mode of quarks. 

Although the relation between the energy scales of con- 
finement and chiral symmetry breaking is interesting, our 
result is not conclusive but suggestive in the present ac- 
curacy. To approach the realistic situation of QCD, we 
would need the reliable chiral extrapolation including the 
dynamical quark effect. 
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